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Abstract 


Results from a study to assess the accuracy of turbulent heating and skin friction prediction 
techniques for hypersonic applications are presented. The study uses the original and a mod- 
ified Baldwin-Lomax turbulence model with a space marching code. Grid converged turbu- 
lent predictions using the wall damping formulation ( original model ) and local damping 
formulation ( modified model ) are compared with experimental data for several flat plates. 
The wall damping and local damping results are similar for hot wall conditions, but differ 
significantly for cold walls, i.e. T w / T r <0.3, with the wall damping heating and skin friction 
10-30% above the local damping results. Furthermore, the local damping predictions have 
reasonable or good agreement with the experimental heating data for all cases. The impact 
of the two formulations on the van Driest damping function and the turbulent eddy viscosity 
distribution for a cold wall case indicate the importance of including temperature gradient 
effects. Grid requirements for accurate turbulent heating predictions are also studied. These 
results indicate that a cell Reynolds number of 1 is required for grid converged heating pre- 
dictions, but coarser grids with a y + less than 2 are adequate for design of hypersonic vehi- 
cles. Based on the results of this study, it is recommended that the local damping formulation 
be used with the Baldwin-Lomax and Cebeci-Smith turbulence models in design and analysis 
ofHyper-X and future hypersonic vehicles. 

1. Introduction 

NASA’s Hyper-X Program was initiated to mature hypersonic air-breathing propulsion, and 
the design and analysis tools required to build an efficient hypersonic cruise vehicle[l]. The 
Hyper-X program will conduct flight tests of an airframe integrated scramjet at Mach num- 
bers of 7 and 10. Successful demonstration of acceleration at Mach 7 and sustained cruise at 
Mach 1 0 will move hypersonic air-breathing vehicle technology from the laboratory into the 
flight environment. The Hyper-X vehicle is designed to have turbulent flow over a significant 
portion of the scramjet engine flowpath. The vehicle structure must be designed to survive 
the increased heating due to this Uirbulent flow, while the overall thrust predictions must con- 
sider the increased skin friction. This connection between the turbulent heating and skin fric- 
tion predictions, and the vehicle design raised concerns regarding the accuracy of turbulent 
flow prediction techniques. To address these concerns, a shidy to assess the accuracy of the 
turbulent heating and skin friction prediction techniques used to design the Hyper-X vehicle 
was undertaken. This study considered both analytical and engineering methods, and Com- 
putational Fluid Dynamics (CFD) codes. The CFD results from the study are presented in 
this report. 

In this study, an upwind, finite volume flow solver, CFL3DE[2], was used to evaluate the 
accuracy of the Baldwin-Lomax algebraic turbulence model [3] for hypersonic applications. 
Shirazi and Truman[4], and Vuong and Coakley[5] have previously compared predictions 
using the Baldwin-Lomax model with experimental data for hypersonic flows. However, the 
wall conditions considered in this study are significantly colder than those from the previous 
studies and span the expected range for the Mach 7 and Mach 10 Hyper-X flights. In the 
present study, grid converged heating and skin friction predictions for flat plates at Mach 
numbers from 6 to 9 were compared with experimental measurements to determine the accu- 
racy of the predictions. In some cases, the CLD solutions significantly overpredicted the tur- 
bulent heating levels. In an effort to improve the agreement, a modification to the algebraic 
turbulence model was incorporated into the CLD code. This modification involved a small 
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change to the model, but had a large impact on the accuracy of the predictions. The modifica- 
tion to the turbulence model and its impact on the heating and skin friction predictions are 
discussed. Also, an assessment of the computational accuracy requirements for turbulent 
heating predictions during vehicle design activities is presented. Some recommendations for 
future design and analysis of hypersonic vehicles are also given. 


2. Turbulence Modeling 

Some of the computational tools used in the design of hypersonic vehicles utilize algebraic 
turbulence models to simulate the effects of turbulence. These algebraic models are easily 
incorporated into CFD codes and require the solution of no additional equations. The Bald- 
win-Lomax model is patterned after the Cebeci- Smith turbulence model[6]. With merged 
inviscid-viscous flowfields, determining the boundary layer edge can be difficult. To avoid 
this difficulty, the Baldwin-Lomax model is normally used in Parabolized Navier- Stokes 
(PNS) and Navier-Stokes (NS) flow solvers. 

The Baldwin-Lomax algebraic turbulence model uses an inner and outer layer formulation. 
In the inner layer, the Uirbulent eddy viscosity, p t , is given by 

(Banner = P(KyD) 2 |“l , (D 


where p is the density, K is the von Karman constant, y is the local distance normal to the 
body surface and Icol is the magnitude of the local vorticity vector. The van Driest damping 
function, D, is given by 


D = 1 


(2) 


where A + is a constant and y + is given by 


y+ = 


Jp^Ky 


Pv 


( 3 ) 


where p w is the density at the wall, p w is the laminar viscosity at the wall and T w is the wall 
shear stress. In the outer layer, the turbulent eddy viscosity is given by 


(Pt>, 


outer 
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( 4 ) 


where K is the Clauser constant and C,.„ is an additional constant. For wall bounded shear 


flows, F wake is given by 


F , = y F 

wake ^liiax n 


( 5 ) 


where y max and F max correspond to the location of the maximum value of the “vorticity” 
function, 

F(y) = y|co|D . (6) 


The Klebanoff intermittency factor is given by 


F Kleb(y) = 1/ 


1 + 5.5 


C Kleby ^| 6 ' 
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where C^i, is a constant. The constants in the model have the following values 

A + =26. 

C cp =1.6 

C Kleb = 0 - 3 

k=0.4 

K=0.0180. 


The Uirbulent eddy viscosity becomes 


ht 


^^irmer’ ^ “ ^cross\ 
^t^outer’ y > ycrossj 


( 8 ) 


where y cross is the smallest value of y at which the inner layer and outer layer formulations 
for the turbulent eddy viscosity are equal. 

The inner layer of the Cebeci-Smith model is similar to eq. [1], but the magnitude of the vor- 
ticity is replaced with the normal gradient of the tangential velocity. In the boundary layer, 
these two quantities are similar. Also, the Cebeci-Smith model uses a different definition for 
y + with 



where p and p are local values of density and laminar viscosity, respectively. The two defini- 
tions of y + , eq.[3] and eq. [9], differ by the factor 



This factor becomes significant in cold wall boundary layers, where the temperature gradient 
at the surface is large. To determine the effect of this factor, the Baldwin-Lomax model in 
CFL3DE was modified to use eq. [9] in place of eq. [3], Using the two definitions for y + as 
guidelines, the van Driest damping function, eq. [2], with eq. [3] for y + is defined as “wall 
damping” and with eq. [9] as “local damping.” An examination of several papers on turbu- 
lence modeling for high speed flows has shown inconsistencies in the manner in which y + is 
defined with these two algebraic turbulence models. For example, Voting and Coakley[5], 
and Stock and Haase[7] use wall damping with both the Cebeci-Smith and Baldwin-Lomax 
models, while Shirazi and Tmman[4] use wall damping for Baldwin-Lomax and local damp- 
ing for Cebeci-Smith. Cheatwood and Thompson[8] use local damping for both turbulence 
models. There appears to be some confusion as to which y + definition should be used in 
these algebraic turbulence models. Therefore, heating and skin friction predictions using 
both definitions for y + are presented in this report. The results will be used to recommend 
guidelines for applying the Baldwin-Lomax turbulence model to hypersonic flows. 


3. Grid Refinement 

Grid refinement was used in this study to evaluate the effect of grid resolution on turbulent 
heating and skin friction predictions, and to obtain grid converged solutions for comparison 
with the experimental data. The grid refinement was perfomted in both the axial and normal 
directions. The nomial grid was created using the grid distribution method of Cheatwood and 
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Thompson[8]. This distribution clusters the grid points to adequately resolve the boundary 
layer gradients without overly resolving the outer shock region. Therefore, grid stretching is 
used near the surface and a uniform spacing is used in the outer region of the grid. At the 
interface of these regions, the grid spacings are forced to match. The fraction of the total 
cells in the stretched region is given by 

F sti = max ^,! - y] , (11) 


where J is the total number of normal cells. The number of stretched cells then becomes 


F su J - 


( 12 ) 


A sinusoidal distribution is used in the recursive formula for the nondimensional height of 
the stretched cells and is given by 


A Yj = 



..sin 


(j - l)rc 
lJsu-i 



i > 


(13) 


where f str is iteratively determined. The total nondimensional distance from the surface to the 
outer edge of a cell is given by 

d4) 

i = 1 

In these calculations, the height of the first cell and the outer boundary are specified. An iter- 
ative procedure is then used to compute the value of f str which yields a nondimensional dis- 
tance from the surface to the outer boundary of 1 . Dimensionality is returned using 

— ^j^standoff’ 0^) 

where y sta ndoff ' s the physical distance between the surface and the outer boundary. This type 
of normal grid distribution is well suited for a grid refinement study involving surface quan- 
tities like heat transfer and skin friction. With relatively small increases in the total number 
of normal grid points, the grid near the surface can be refined as needed, while an adequate 
uniform spacing in the outer region remains relatively unchanged. 

The axial grid starts slightly downstream of the leading edge of the flat plate, with the axial 
grid step growing by a factor of 1.03 until a specified maximum axial step is reached. During 
the grid refinement, the first axial plane is moved forward and the maximum axial step is 
reduced. 


Cheatwood and Thompson[8], and Weilmuenster and Gnoffo[9] recommend a cell Reynolds 
number of 1 for turbulent calculations, with the cell Reynolds number defined as 


Re 


cell 


Pt a t A y i 

hi 


(16) 


where p a 1 and P| are the density, speed of sound and laminar viscosity, respectively, at the 
first cell and Ay! is the cell height. However, the standard measure of grid resolution for vis- 
cous flows uses the value of y + at the first point off the surface. To accurately resolve the 
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temperature profile, Fletcher[10] recommends a y + less than 2. The grid refinement study 
will address the required cell Reynolds number and y + for efficient and accurate turbulent 
heating and skin friction predictions. 

4. Computational Procedure 

Previous studies using this upwind, finite volume flow solver have demonstrated its ability to 
accurately predict the flowfield on sharp cones[l 1], sharp nose hypersonic forebodies [12] 
and hypersonic forebody/inlet flowpaths[13,14]. Results with grid refinement have shown 
reasonable agreement with experimental heat transfer measurements on forebodies at Mach 
numbers from 11 to 19 [13,14], 

The PNS solution was started using conical flow at the first axial station. This required that 
the normal grid points in the first cell extrapolate to the leading edge of the flat plate, with the 
first cell located slightly downstream of the leading edge, as shown in Figure 1 . After a coni- 
cal solution was obtained, the marching solution was initiated. The grid was constmcted with 
the outer boundary located to capture the Mach line originating from the nose of the flat 
plate. 

Transition to turbulence was calculated using the exponential stream wise intermittency func- 
tion of Dhawan and Narasimha[15], with the onset and end of transition based upon the 
experimental data. Turbulent flow was modeled using the original and modified Baldwin- 
Lomax model as discussed previously. 

Solutions were considered to be grid converged when the heat transfer and skin friction dis- 
tributions were unchanged with further grid refinement. Grid converged solutions were first 
obtained using local damping and then, wall damping solutions were obtained on the finest 
grid. Results from a boundary layer code [16] are also presented for reference. The boundary 
layer code uses the Cebeci-Smith turbulence model with local damping in the van Driest 
damping function. 
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Figure 1. Conical Grid at Nose of Flat Plate 


5. Results 

5.1 Experimental Database 

A large database of experimental heat transfer and skin friction measurements on simple 
configurations is available for assessing the accuracy of hypersonic turbulent flow prediction 
techniques. The data covers a wide range of wall conditions. The expected wall temperatures 
on the Mach 7 and Mach 10 flight vehicles vary with the tungsten forebody as low as 1000° 
R and the themial protection system (TPS) tiles as high as 2500° R. Relative to the flight 
total temperatures, these wall temperatures are rather cold. For purposes of this study, a cold 
wall is defined as having a T w / T t less than 0.3, where T w is the wall temperature and T t is 
the freestream total temperature. The majority of the results presented in this report are for 
cold walls. The experimental datasets are from Bertram, Cary and Whitehead [17]; Bertram 
and Neal[18]; and Wallace[19], The flow conditions for all the cases are given in Table 1. 
The Bertram, Cary and Whitehead, and Bertram and Neal datasets are from Mach 6 flat plate 
experiments conducted in the 20-Inch hypersonic tunnel at NASA Langley Research Center. 
The data consist of Stanton number distributions. One of the cases has a cold wall, while the 
remaining cases are near adiabatic wall conditions. Also, one of the cases has the flat plate 
oriented at an angle of attack of 15°. The Wallace flat plate data covers a Mach number range 
from 6.4 to 8.8. All of the cases have a cold wall, with a total enthalpy ratio, H w / H t , varying 
from 0. 1 1 to 0.26. With the perfect gas assumption from the experiments, the total enthalpy 
ratio, H w / H t , is equivalent to the temperature ratio, T w / T t . The data consist of Stanton 
number and skin friction coefficient distributions. The three Wallace cases considered in this 
study provide a good sampling of the range of Mach number and wall temperature from the 
Wallace dataset. The range of wall conditions from these experimental datasets span the 
expected wall conditions from the Hyper-X flights (0.1 1 < T w / T t < 0.57). 
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Table 1. Flow Conditions for Flat Plate Cases 


Reference 

Case 

Description 

M oo 

Too (K) 

R e oo (l /m ) 

T w (K) 

Data 

Bertram, Cary and 
Whitehead[ 17] 

T w / T t = 0.2 

6.0 

65.04 


106.67 

N st 

T w / T t = 0.6 

6.0 

65.04 


320.0 

N st 

Bertram and 
Neal[18] 

T w / T t = 0.65 
a = 0.0° 

6.0 

63.01 

26.38xl0 6 

335.84 

N st 

T w / T t = 0.65 
a = 15.0° 

6.0 

63.01 

26.38xl0 6 

335.84 

N st 

Wallace[19] 

H w /H t = 0.11 

7.12 

240.0 


299.0 

N st 

H w /H t = 0.19 

8.80 

93.99 


296.3 

N s t> Cf 

H w /H t = 0.26 

7.00 

105.0 


296.0 

N st , C f 


5.2 Mach 6 Cases 

The grid dimensions and spacings, and transition onset and end locations for the Mach 6 
cases are given in Table 2. The grid converged solution is with grid 1 and the finest grid con- 
sidered is grid 2. The Re ce jj and y + for each grid are also given. In some instances, the maxi- 
mum axial grid spacing, Ax max , was not changed when the normal grid was refined. Only 
heat transfer measurements were reported for the Mach 6 flat plate tests with the experimen- 
tal data consisting of Stanton number distributions, where the measured heat transfer has 
been normalized using the edge conditions. For the angle of attack case, the edge conditions 
are the post-shock conditions, i.e. M edge =4.0, and not the freestream conditions. 


Table 2. Mach 6 Cases 


Case 

Description 

Grid 

R e cell 

+ 

yi 

Grid Dimensions 
(Normal x Axial) 

Ay x 

(in) 

Ax max 

(in) 

Transition 
onset - end 

T w /T t = 0.20 

1 

2.0 

.15 

161 x 191 

.38x10"' 


4.8" - 10.4" 


2 

1.0 

.07 

201x311 

.19x10 4 


T w / T t = 0.60 

1 


.07 

121 x 191 

Hill 

.10 

3.75" -7.5" 


2 


.03 

161 x 191 


.10 

T w / T t = 0.65 
a = 0.0° 

1 

2 

1.0 

0.4 

.07 

.03 

121 x 191 
161 x 191 

m 

.10 

.10 

3.75" -7.5" 

T w / T t = 0.65 

1 

2.5 

.13 

161 x 191 



3.15" -5.25" 

a = 15.0° 

2 

1.0 

.05 

201 x 301 




The results for the cold wall case, T w /T t =0.2, are plotted using the running Reynolds number, 
R x , in Figures 2 and 3. Comparisons of CFD and boundary layer heating predictions with 
experimental data are shown in Figure 2. The skin friction results are shown in Figure 3. The 
local damping results with grid 2, the finest grid, differ by less than 0.3% from the grid 1 
results, indicating the solution is grid converged. The grid 1 and grid 2 solutions have 95 and 
120 points in the boundary layer, respectively. The local damping predictions are in reason- 
able agreement with the data with a Root Mean Square (RMS) deviation of 15% from the 
turbulent data. The wall damping predictions are about 23% above the local damping results 
and have an RMS deviation of 41% from the turbulent data. The boundary layer heating and 
skin friction results are in excellent agreement with the local damping predictions. The 
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results for the hot wall case, T w /T t =0.6, are shown in Figures 4 and 5. Comparisons of CFD 
and boundary layer heating predictions with experimental data are shown in Figure 4. The 
skin friction results are shown in Figure 5. The CFD predictions using local and wall damp- 
ing differ by only 2% and have good agreement with the data with an RMS deviation of 8% 
from the turbulent data. The boundary layer results are 4% higher than the CFD. The grid 1 
and grid 2 solutions have 65 and 94 points in the boundary layer, respectively. These two 
cases show the effect of wall temperature with the local damping and wall damping formula- 
tions of the Baldwin-Lomax turbulence model. For hot walls with small temperature gradi- 
ents, die two formulations give slightly different results; but for cold walls with large 
temperature gradients, the results differ significantly. Only the local damping results have 
good or reasonable agreement with the data for both cases. 

Comparisons of CFD and boundary layer heating predictions with experimental data for the 
other Mach 6 cases are shown in Figures 6 and 7. The 0.0° angle of attack CFD results in 
Figure 6 have fair agreement with the data with an RMS deviation of 18% with local damp- 
ing and 21% with wall damping. The 15.0° angle of attack CFD results in Figure 7 are plot- 
ted using local edge conditions and have good agreement with the data with an RMS 
deviation of 7% with both local and wall damping. The good agreement at 15.0° angle of 
attack and the large overshoot in the 0.0° angle of attack data at the end of transition suggest 
a possible tunnel-induced disturbance causing the less than good agreement in Figure 6. For 
both cases, the boundary layer results were 5% above the CFD predictions. These hot wall 
CFD results are similar to the previous hot wall case with the two formulations yielding only 
slight differences. 



Figure 2. Mach 6 Cold Wall Heating Comparisons 
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Figure 3. Mach 6 Cold Wail Skin Friction Comparisons 



°o,X 


Figure 4. Mach 6 Hot Wall Heating Comparisons 
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Figure 5. Mach 6 Hot Wall Skin Friction Comparisons 



Figure 6. Mach 6 Hot Wall Heating Comparisons 
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Figure 7. Mach 6 Hot Wall, Angle of Attack Heating Comparisons 


5.3 Wallace Cold Wall Cases 

The Mach 6 results indicate the turbulent heating predictions with local damping and wall 
damping differ significantly under cold wall conditions. Because accurate cold wall heating 
predictions are crucial to the themial analysis of the Hyper-X flight vehicles and the related 
TPS sizing, cold wall effects were further studied with several cases from the Wallace 
dataset[19]. The grid dimensions and spacings, and transition onset and end locations for the 
Wallace cases are given in Table 3. The Re cc u and y + for each grid are also given. The exper- 
imental data consist of Stanton number and skin friction coefficient distributions, where the 
measured quantities have been nomialized using the edge conditions. 


Table 3. Wallace Cold Wall Cases 


Case 

Description 

Grid 

R e cell 

+ 

yt 

Grid Dimensions 
(Normal x Axial) 

(in) 

Ax max 

(in) 

Transition 
onset - end 

Mach 7.12 

1 

2.2 

.18 

161x211 

.5x10 4 

mm 

5.8"- 11.5" 

H w / H t = 0. 1 1 

2 

1.0 

.08 

201 x 321 

.23x10 4 


Mach 8.80 

1 

2.2 

.17 

161x211 

ISjH 

wa 

3.9" - 8.0" 

H w / H t = 0.19 

2 

1.0 

.08 

201 x 321 



Mach 7.00 

1 

2.3 


161x211 

.3x10 4 

■a 

2.0" -4.5" 

H w / H t = 0.26 

2 

1.0 


201 x 321 

.14x10 4 

m 


11 





































The results for the Wallace cases are plotted against the running Reynolds number, R x , in 
Figures 8-13. For each case, the grid 1 and grid 2 solutions have 96 and 126 points in the 
boundary layer, respectively. Comparisons of CFD and boundary layer heating predictions 
with the Mach 7.12 experimental data are shown in Figure 8. The skin friction results for that 
case are shown in Figure 9. The local damping predictions have good agreement with the 
heating data with an RMS deviation of 5%. The wall damping heating and skin friction pre- 
dictions are 32% above the local damping results and have an RMS deviation of 39% with 
the heating data. Comparisons of CFD and boundary layer predictions with the Mach 8.80 
experimental data are shown in Figures 10 and 11. The local damping predictions have rea- 
sonable agreement with the heating data and fair agreement with the skin friction data, with 
RMS deviations of 1 1% and 19%, respectively. The wall damping predictions are 19% above 
the local damping results and have RMS deviations of 31% with the heating data and 40% 
with the skin friction data. Comparisons of CFD and boundary layer predictions with the 
Mach 7.0 experimental data are shown in Figures 12 and 13. The local damping predictions 
have reasonable agreement with the data with RMS deviations of 14% with the heating data 
and 12% with the friction data. The wall damping predictions are 12% above the local damp- 
ing results and have RMS deviations of 28% with the heating data and 22% with the friction 
data. For all cases, the boundary layer results and the local damping predictions are within 
3%. 

Profiles from the Mach 7.12 solutions at the last axial station have been studied to determine 
the impact of the two damping formulations on the flowfield properties and the turbulent 
eddy viscosity distribution. The full temperature profiles are shown in Figure 14, with a 
blow-up of the region near the surface in Figure 15. The agreement between the local damp- 
ing solution and the boundary layer results is good through the inner layer region, but they 
differ somewhat in the outer portion of the boundary layer. This may be attributable to the 
differences in the modeling of the outer layer region between the Baldwin-Lomax and 
Cebeci-Smith models. Also, the wall damping solution has a thicker boundary layer than the 
local damping solution. Figures 14 and 15 have large temperature gradients near the surface, 
with a peak in temperature at approximately 0.005 inches from the surface. The impact of 
this temperature gradient on the van Driest damping function, eq. [2], with local damping 
and wall damping is shown in Figure 16. The damping function with either formulation 
becomes 1.0 above a height of 0.06 inches, but the damping function's exponential temi 
decays significantly slower with local damping. The turbulent eddy viscosity profiles near 
the surface are shown in Figure 17. A slower initial growth in p, is present with local damp- 
ing due to the increased damping in the first 0.06 inches shown in Figure 16. The turbulent 
eddy viscosity profiles in the inner and outer layers are shown in Figure 18. The wall damp- 
ing formulation has larger values of eddy viscosity than the local damping, with the peak val- 
ues in eddy viscosity 37% higher with wall damping. The crossover point between the inner 
and outer layer formulations occurs at about the same location with local damping and wall 
damping, but the wall damping eddy viscosity is already 33% above the local damping value 
at that point. 

These results show a growing difference between the CFD predictions with local damping 
and wall damping as T w /T t or H w /H t is decreased. In all cases, the local damping predictions 
had good or reasonable agreement with the experimental heating data, while the wall damp- 
ing predictions were at least 25% above the data. In the lowest T w /T t case, the wall damping 
predictions were 39% higher than the data. The local damping predictions had reasonable or 
fair agreement with the experimental skin friction data, and again the wall damping predic- 
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tions were at least 22% above the data. A 30% overprediction of turbulent heating and sk in 
friction will adversely affect the TPS sizing and drag estimates for the Hyper-X flight vehi- 
cles. These results show the importance of temperature gradient effects in turbulence model- 
ing for hypersonic flows. 



Figure 8. Mach 7.12 Heating Comparisons 
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Figure 9. Mach 7.12 Skin Friction Comparisons 



Figure 10. Mach 8.80 Heating Comparisons 
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Figure 13. Mach 7.0 Skin Friction Comparisons 
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Figure 15. Mach 7.12 Temperature Profiles Near the Surface 








5.4 Computational Requirements for Accurate Heating Predictions 

The cold wall cases required a Re cell of 2 for grid converged solutions, while the hot wall 
cases required a Re cel j of 1 . However, these computational grids are not practical for analyz- 
ing large scale hypersonic vehicles because of the fine normal and axial grid spacings. To 
assess the accuracy requirements for CFD turbulent heating predictions, coarser grids were 
used to compute the flowfields for the Wallace Mach 7.12 and the Bertram, Cary and White- 
head T w /T t =0.6 cases. The grid dimensions and spacings, and Re ccl ], y + and number of points 
in the boundary layer are given in Table 4, with the grids numbered from the coarsest to the 
finest. The percentage change from the turbulent heating levels with the finest grid are given 
in the last column. In both cases, the coarsest grid has less than 20% of the total grid points 
of the finest grid and 30% of the grid points in the boundary layer. The results for the Wallace 
case are shown in Figures 19-22. The heating predictions in Figure 19 show a 5.8% differ- 
ence with grid 1 and a 2.8% difference with grid 2. The skin friction results in Figure 20 
show a smaller change, with the grid 1 results 1.8% below the finest grid. The temperature 
profiles at the last axial station in Figure 21 have excellent agreement. The symbols at each 
cell center also show the refinement of the grid was concentrated near the surface, while the 
grid near the boundary layer edge had much less refinement. The temperature profiles in Fig- 
ure 22 show small differences at the surface due to the coarsening. The y + of the coarsest 
grid satisfies Fletcher’s recommendation! 10] of y + less than 2. The results for the Bertram, 
Cary and Whitehead case are shown in Figures 23-26. The heating predictions in Figure 23 
have a 5.7% difference with grid 1 and a 2.9% difference with grid 2. The skin friction 
results in Figure 24 again have a smaller change, with grid 1 results differing by 0.7% from 
the finest grid. As in the previous case, the temperature profiles at the last axial station in Fig- 
ures 25 and 26 have excellent agreement. In this case, the y + of the coarsest grid is smaller 
than with the cold wall. These results indicate that while a Re ce jp 1 .0 may be required for 
grid converged heating predictions for these flat plate cases, this requirement is overly 
restrictive for design activities. For design purposes, a 6% underprediction in heating levels 
is acceptable and the reduction in computational time with the smaller grids enhances the 
overall efficiency of the analysis. Because of the limited nature of this study, this conclusion 
is limited to sharp nose geometries, and blunt geometries with a localized stagnation point 
may require a Re ce ||= 1 .0. 


Table 4. Heating Prediction Accuracy Requirement Cases 


Case Description 

Grid 

Re cell 

+ 

yt 

# pts 

Grid 

Dimensions 

A Yl 

(in) 

Ax max 

(in) 

A N St ~ 

Wallace 

1 

17.2 

1.46 

37 

81 x91 

.50x10 3 

.502 

5.8% 

Mach 7.12 

2 

9.4 

0.77 

52 

101 x 121 

.24x10 3 

.335 

2.8% 

H w / H t = 0. 1 1 

3 

1.0 

0.08 

126 

201 x 321 

.23x10 4 

.100 

- 

Bertram, Cary and 
Whitehead 
Mach 6 
T w /T t = 0.6 

1 

7.1 


35 

61 x81 

.60x10 3 

.3455 

5.7% 

2 

3.6 


51 

81 x 101 

.30x10 3 

.2255 

2.9% 

3 

0.5 


120 

161 x 191 

.38x10 4 

.1000 
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Figure 19. Mach 7.12 Heating Comparisons with Coarser Grids 




Figure 20. Mach 7.12 Skin Friction Comparisons with Coarser Grids 
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Figure 21. Mach 7.12 Temperature Profiles with Coarser Grids 
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Figure 22. Mach 7.12 Temperature Profiles Near the Surface 
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Figure 24. Mach 6 Skin Friction Comparisons with Coarser Grids 
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6. Conclusions 


Results from a study to assess the accuracy of CFD turbulent heating and skin friction pre- 
diction techniques for hypersonic applications have been presented. The study used the orig- 
inal and a modified Baldwin-Lomax turbulence model with a PNS code. Grid converged 
turbulent predictions using the wall damping formulation (original model) and local damping 
formulation (modified model) were compared with experimental data for several flat plates. 
The wall damping and local damping results were similar for hot wall conditions, but dif- 
fered significantly for cold walls with the wall damping heating and skin friction 10-30% 
above the local damping values. Furthermore, the local damping predictions had reasonable 
or good agreement with the experimental heating data and reasonable or fair agreement with 
the experimental skin friction data. The impact of the two damping formulations on the van 
Driest damping function and the turbulent eddy viscosity distribution was shown for a cold 
wall case, indicating the importance of temperature gradient effects. Grid requirements for 
accurate turbulent heating predictions were also studied by using coarser grids to obtain pre- 
dictions for a cold wall and a hot wall case. These results indicated that a cell Reynolds num- 
ber of 1 may be required for grid converged solutions, but coarser grids are adequate for 
design purposes. The cold wall results supported the less restrictive requirement of y + less 
than 2. 

The comparisons between CFD predictions and experimental data presented in this report 
are for hypersonic flat plate flowfields. In all cases, the local damping formulation had as 
good or better agreement with the data than the wall damping. At the cold wall conditions, 
the wall damping overpredicted the turbulent heating data by as much as 40%. Therefore, it 
is recommended that the local damping formulation be used with the Baldwin-Lomax and 
Cebeci-Smith turbulence models for design and analysis of the Hyper-X flight vehicles. 
Because the Hyper-X design uses forebody compression surfaces on the engine flowpath, 
further study is required to determine the impact of local damping and wall damping in 
regions of isentropic and discrete compression. But given the current results, the use of local 
damping with algebraic turbulence models is recommended for obtaining accurate turbulent 
heating and skin friction predictions. 
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